Enhancing the spin-transfer torque through proximity of quantum well states 



loannis Theodonis, 1, 2, Q Alan Kalitsov, 1 and Nicholas Kioussis 1 

1 Department of Physics, California State University, Northridge, CA 91330-8268 
A Department of Physics, National Technical University, GR-15773, Zografou, Athens, Greece 

(Dated: February 6, 2008) 

We predict that the spin-transfer, T t n, and field-like, Ti,x, components of the local spin torque are 
dramatically enhanced in double-barrier magnetic tunnel junctions. The spin-mixing enhancement 
is due to the energetic proximity of majority and minority quantum well states (QWS) of different 
quantum numbers within the bias window. T$ n exhibits a switch-on and switch-off step-like bias 
behavior when spin polarized QWS enter the bias window or exit the energy band, while T it ±, 
changes sign between switch-on biases. The net T± exhibits an anomalous angular behavior due to 
the bias interplay of the bilinear and biquadratic effective exchange couplings. 
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I. INTRODUCTION 

Spintronics involve the exploitation of the quantum- 
mechanical spin degree of freedom to provide new func- 
tionalities beyond conventional electronics [H . One effect 
that has its roots on the electron's spin, is the torque 
exerted on the magnetization of a nanometer-scale free 
ferromagnet (FM) by a spin-polarized current, originat- 
ing from a preceding non-collinear pinned FM0, y, 0j- 
This torque can be decomposed into a field-like and a 
spin-transfer component both orthogonal to the mag- 
netic moment of the free FM, but with different influ- 
ence on its dynamics [6|. Recent ferromagnetic resonance 
experiments, provide a useful tool to study the role of 
each component 0, Q. At sufficiently high current den- 
sities, the spin-transfer torque leads to current-induced 
magnetization switching (CIMS) [1, EH EH • Reduction of 
the high critical current for CIMS is necessary for spin- 
transfer controlled magnetic memories fl2l]. 

Double-barrier magnetic tunnel junctions (DBMTJ) 
consist of a central metallic layer between two insu- 
lating barriers and two FM electrodes. The tunneling 
magnetoresistance (TMR) can be dramatically enhanced 
in collinear DBMTJ by the presence of quantum well 
states (QWS) under appropriate resonant conditions [l3l 
[LH HH. Recently, the discrete energy spectrum of 
FM nanoparticles has been shown to enhance spin ac- 
cumulation and to control the bias dependence of the 
TMR[i3|. While the physics in collinear DBMTJ has 
been studied extensively [H, EH EE El]) the effect of 
spin-polarized QWS (SPQWS) on the spin-torque in non- 
collinear DBMTJ remains an unexplored area thus far. 

The objective of this work is to present for the first time 
a study of the effect of SPQWS on both the spin-transfer 
T|| and field-like T± components of the spin torque under 
external bias. The calculations are based on the tight- 
binding method and the non-equilibrium Keldysh formal- 
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FIG. 1: (Color Online) Schematic of the DBMTJ consisting 
of a FM central wire of Nc atomic sites, connected to left and 
right FM leads through the tunneling barriers I of N sites. 
The spin-transfer (field-like), T t ii (T;,j_) components of the 
torque lies in the -x (y) directions. 



ism. We predict that both components of the local spin 
torque can be dramatically enhanced when majority and 
minority QWS energies of different quantum numbers are 
in close proximity and lie within the bias window. It 
should be emphasized that the local-spin-torque enhance- 
ment is not associated with an enhancement of the corre- 
sponding spin-polarized currents. The low-temperature 
bias dependence of the local spin-transfer torque, ex- 
hibits a switch-on and switch-off step- like behavior when 
the SPQWS enter the bias window or exit the energy 
band, respectively, similar to that of the spin-polarized 
currents. On the other hand, T^± changes sign between 
the majority and minority switch-on values of bias. We 
demonstrate that the bias behavior of Tj n and T^± can 
be derived analytically using a single-site central FM re- 
gion. The net Tj_, pertinent to the non-equilibrium in- 
ter layer exchange coupling, Exc, exhibits an anomalous 
angular behavior due to the interplay of the bilinear and 
biquadratic effective exchange couplings which have dif- 
ferent bias behavior. 

The structure of this paper is as follows. In Sec. II, 
we introduce the basic model and outline the computa- 
tional approach of the non-equilibrium spin torque. In 
Sec. Ill, we present and discuss the results of the calcu- 
lations. Finally, Sec. IV includes a brief statement of the 
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II. METHOD - FORMALISM 

FigCD shows the one-dimensional FM/I/FM/I/FM 
DBMTJ system, consisting of a central (C) FM nanowire 
containing Nc atomic sites (AS) connected to the left (L) 
and right (R) FM electrodes through two thin symmet- 
ric non-magnetic tunneling barrier nanowires / of N = 2 
AS. The magnetization of the central FM, Mp, is along 
the z axis of the coordinate system shown in Fig. [1] The 
magnetization of the FM leads M^^) lies in the x — z 
plane, i.e. it is rotated by angle around the wire axis 
y. The Hamiltonian of the system is 

H = H L + H M + H R + H L , M + H RM + hc, (1) 

where Hm = Hj + He + Hj c + he is the Hamil- 
tonian of the middle (M) multi-layer //FM// region, 
and Hl/rym> Hi.c are the coupling Hamiltonians at the 
L(R)/M and I/C interfaces, respectively. The Hamilto- 
nian for each FM region H a , a= L, R, and C, is described 
by a one-dimensional single-orbital tight-binding model 
neglecting the in-plane k-dependence, which includes a 
nearest-neighbor (NN) spin-independent hopping term, 
t a , and a spin-dependent on-site energy term, e^, i.e., 



t a c\ci. 



he. 



(2) 



The Hamiltonian, Hj, for the barriers is identical to 
H a , but where one replaces the hopping term t a with tj 
and the spin-dependent on-site energy e£ with the spin- 
independent £/. The coupling Hamiltonian of the FM/I 
interface is H a j — t a /ic\.ci. The exchange-splitting, 
A Q = 4 — 4j is identical in all FMs, Ai = 0, 4 — E F = 
0.318 eV, 4 " e f = 0.736 eV, and e x - E F = 6.5 eV, 
where Ep is the Fermi energy. The t a — 0.4 eV in all 
FMs, ti = 1 eV in both insulators, and t a /j = 0.5 eV in 
all FM/I interfaces, consistent with the ab initio values 
for ID Co FM nanowires [Hj]. 

We extend Datta's formalism [l9[ to the case of non- 
collincar systems, where the scalar Green functions are 
replaced with 2x2 matrices in spin space. For this pur- 
pose, H a can be expressed in the form, H a = H a + 5H a , 
where 
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where / is a 2x2 unit matrix. Following Datta [19| we find 
that the retarded Green function of the middle region M 
is 
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AI 



(6) 



where E is the one-electron electron energy, Hm and 
are the (2Nm x 2Nm) Hamiltonian and self-energy 
matrices, respectively, and Nm — 2N+Nc is the number 
of atomic sites in the middle region. The only non-zero 
elements of ^l(r) are the (2 x 2) self-energy matrices at 
the interfacial sites 



X L (R)(0)=t 2 a/I g L (R)(0), 



(7) 



where gi,(R){0) are the retarded surface 2x2 Green's 
function matrices of the isolated L(R) lead, determined 
from Eq. O The non-equilibrium Green's functions can 
be determined by solving the kinetic equation [T|| 



G M — iGu^ < G M , 



(8) 



where S< = f L (t[ - t L ) + f R (t R - t R ), is the non- 
equilibrium self-energy matrix and /l(r) are the Fermi- 
Dirac distribution functions of the L(R) leads. 

The local spin-transfer torque exerted on the local 
moment at site i in the central FM region is(20j 



where 



T< = —V • I 



tc 



M _ r(«) 
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dE, 



(9) 



(10) 



is the spin current between NN sites [2fJ, and cr = 
(a~ x ,(Ty,cr z ) is a vector of the Pauli matrices. The 
field-like, Ti t ±, and spin-transfer, T^m, components of 
the local spin torque, shown in Fig. [TJ are along the 
M c x (Ml(r) x M c ) and M c x Ml(r) directions, respec- 
tively. Here, Mc and M L ^ are unit vectors along the 
magnetization of the free C and pinned L(R) FM regions, 
respectively. 



describes the spin-average part of H a and 

«f Q = §(4-4) 



(4) 



is the spin-splitting part of H a . The one-electron 
Schrodinger equation for the retarded Green function, 
g°'° , in each isolated semi-infinite ferromagnetic lead be- 



lli. RESULTS - DISCUSSION 

The majority-(full triangles) and minority- (open tri- 
angles) QWS energies relative to the Fermi energy, E n , 
as a function of the thickness, Nc, of the central FM 
wire are shown in Fig. [2] for QWS between -0.5 eV and 
0.5 eV. The numbers next to each series of data points, 
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FIG. 2: (Color online) SPQWS energy positions E" as a func- 
tion of the number of atomic sites, N c , in the central FM re- 
gion for zero bias and 9 — 0. The bottom of the majority 
(minority) conduction band of the leads is denoted by -Ejj 
and A is the exchange splitting, denoted by the dashed hori- 
zontal arrow. At finite bias V , the chemical potentials of the 
L,R leads are shifted by eV = /j.r — fiL around the Fermi 
energy. 



indicate the quantum number, n a = V , 2°", . . . , N% , of 
the SPQWS. The dashed curves denote the SPQWS en- 
ergies, Ef = e a + 2tcos{n' J ir/(Nc + 1)) of the isolated 
central FM wire. The coupling of the central region to 
the FM leads results in a shift and a broadening of the 
SPQWS energies. The bottom of the majority (minority) 
conduction band of the leads, at zero bias, is indicated 
by £]P. Note, that for N c = 7AS the n T = 3 T majority- 
and = l- 1 minority QWS are in very close proximity 
and they are very close to the Fermi energy. Under ap- 
plied bias V, only the SPQWS with energies E n " that lie 
within the bias window from /zl = — ^g- to fiji = +^f, 
denoted by the shaded area in Fig. [21 contribute to the 
resonant tunneling. 

In Figs. [3^ and[3j3 we show the perpendicular, Ti t ±, 
and parallel, T^u, components of the local spin torque 
on the first site (i = 1) in the central FM region next 
to the left FM/I interface, as a function of the thickness 
N c of the central FM region. Both local spin torque 
components are calculated at T=5K and V=0.1V, for 
the almost parallel 8 — 7r/10, perpendicular 9 = ir/2 and 
almost antiparallel 9 = 97r/10 configurations. We find 
that Ti t ± and Tj n are strongly enhanced for N c = 7 AS 
by about one and two orders of magnitude, respectively. 
The local spin-torque enhancement persists even for very 
small angular deviations from the parallel(circles) and 
antiparallel(triangles) configurations. Interestingly, T^j_ 
for N c = 7 changes sign with increasing 9, in contrast to 



T, 




FIG. 3: (Color online) Low-temperature (T=5 K) and low- 
bias (V = 0.1 V) (a) perpendicular, Ti,x, and (b) parallel, 
Ti M, components of the local spin torque versus the number of 
atomic sites, N c , of the central FM region. Both components 
of the local spin torques are calculated on the first site in the 
central FM region next to the left FM/I interface for three 
values of the angle 8 of n/10, tt/2, and 97r/10, respectively. 
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FIG. 4: (Color online) Angular dependence of T it ± (black 
squares) and T it \\ (red circles) for N c — 4AS in panels (a) and 
(b), and for N c = 7 AS in panels (c) and (d), respectively. The 
local spin torque is evaluated on the first site in the central 
FM region next to the left FM/I interface at T = 5K and 
V = 0.1V. 



pcndcnce of the perpendicular, T^±, (black squares) and 
parallel Tj i (red circles) components of the local spin 
torque for N c = 4 AS and N c = 7AS, respectively. The 
spin torques are evaluated on the first site in the central 
FM region next to the left FM/I interface at T=5K and 
V=0.1V. In both cases (Figs. \kp andgji) Z^, exhibits a 



In Figs. S^IUi and [3Jd|3J1 we display the angular de- sinusoidal angular dependence similar to that in a single 
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FIG. 5: (Color online) Perpendicular, T it ± (black squares, 
and parallel, J^n (red circles), components of spin torque as 
a function of site i in the central FM region, for N c = 4AS 
in panels (a) and (b) and for N c = 7AS in panels (c),(d), 
respectively. The local components are calculated at T=5K, 
0.1 V bias, and 6 = tt/2. 



MTJ0. On the other hand, when the enhancement con- 
ditions are fulfilled for N c = 7 AS (Fig. gfc), T ly± exhibits 
a non-sinusoidal angular dependence changing sign in the 
[0, 7r] interval. 

In Figs. 13 a), (c) and[5jb),(d) we display the perpen- 
dicular, Ti t x, (black squares) and parallel Tj n (red cir- 
cles) components of the local spin torque as a function of 
site i in the central FM, for N c = 4AS and N c = 7AS, 
respectively. The spin torques are evaluated at T=5K, 
V=0.1V, and 6 = tt/2. Interestingly both T i<± and T it \\ 
oscillate around zero as function of atomic site i due to 
the electron precession in the central FM. The number 
of nodes increases as the width of the FM quantum well 



increases. The enhancement of T^x and Tj i for iV c = 7 
AS holds for all sites i. 

In Figs. [6^i and [6Jd, we display the low-temperature 
(T = 5K) bias dependence of the spin-polarized cur- 
rents, P and I 1 , for N c = 4AS and N c = 7AS, re- 
spectively, and for 9 = |. The spin-polarized cur- 
rents switch on at V° n = 2\E n " — E F \, when the n a 
QWS enters the bias energy window. For N c = 7 AS, 
the switch on for both spin-polarized currents occur at 
about the same bias, Vj n « V} n , due to the fact that 
\E lL — Ep\ k \E^ - E F \ in Fig. [2J Both currents de- 
crease with increasing bias because the density of states 
of the minority band in the leads at E n decreases. At 



the critical bias, V l ff = 2(E n + 

gies, E l1 and i? 3 ' are shifted below the bottom of the 
minority band of the lead, and hence the minority con- 
tribution to the spin-polarized currents is switched off. 
In Figs. [5J: and [SJi we display the low-temperature 
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FIG. 6: (Color online) Low-temperature (T=5K) bias de- 
pendence of the majority (green symbols) and minority (blue 
symbols) currents in the central FM for (a) N c — 4 AS and 
(b) N c = 7AS, respectively, and for 9 — tt/2. Bias depen- 
dence of T^n (black symbols) and Ti t ± (red symbols) on the 
first site in the central FM, for (c) N c = 4 AS and (d) N c = 7. 
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denote the switch-on (-off) bias. 



(T = 5K) bias behavior of the parallel, Tf \\, (black 
squares) and perpendicular T^±_ (red circles) components 
of the local spin torque on the first site in the central FM, 
for = | and for N c = 4 AS and N c — 7 AS, respectively. 
The local spin-transfer component, T^u, exhibits a switch 
on bias behavior at Vj n and at V} n , similar to that of the 
spin polarized currents in[6K. On the other hand, Ti t ±, 
which is non-zero for zero bias, displays a non-monotonic 
bias dependence, changing sign between Vj n and V} n , 
similar to that of the exchange field in quantum dots 
connected to FM leads [13] • It is important to note that 
both Tj ii and T^x are strongly enhanced for iV c = 7 AS 
in Fig. [Hji, even though the corresponding spin-polarized 
currents are smaller than those for N c = 4 AS. Thus, the 
enhancement of the local spin torque is not associated 
with a corresponding enhancement of the spin-polarized 
currents. This result clearly demonstrates that the un- 
derlying mechanism that controls the local spin-transfer 
torque enhancement is the close proximity of the majority 
and minority QWS energies of different quantum num- 
ber, E n sa E n , within the bias energy window. This in 
turn enhances the spin mixing a <-» a in the central FM, 
when electrons tunnel resonantly through the SPQWS. 
The enhancement of the local spin-transfer torque is in- 
dependent of the parity of the QWS wavefunctions. This 
spin-transfer torque enhancement may have technologi- 
cal applications since the CIMS may be facilitated under 
such conditions. 

In order to elucidate the role of the SPQWS on the en- 
hancement and bias dependence of the local spin torque, 
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FIG. 7: (Color Online) (al),(a2) Energy dependence of the 
real (red squares) and imaginary (black squares) parts of 
G <n for (£ T << E } ) and (£ T w E v ), respectively;(bl) and 
(b2) bias dependence of T^n (blue squares) and of T iy ± (green 
squares) corresponding to (al) and (a2), respectively. 



we have used a simple model of a central FM region 
consisting of a single-site coupled weakly to the leads. 
We have shown that to leading order in the coupling pa- 
rameter, tc, the on-site spin off-diagonal non-equilibrium 
Greens function is 



a=LM 



(11) 



Here, E£l = t^Tr[Nl — N^]sin(9) is the self-energy due 

1(1) 

to the leads, Na is the surface density of states of 
the L(R) leads, which for simplicity are taken to be en- 
ergy independent. In the weak coupling regime, the re- 
tarded (advanced) Greens functions of the coupled sys- 
tem, G"? a \ , can be be approximated with those of the un- 

coupled system #f a) , i.e. G% a) w g% a) = [E-E°±irj\- X , 
where rj is taken to be spin-independent. Therefore, the 
spin-transfer torque components, determined by the non- 
equilibrium on-site magnetization [2l| . are 



Ti t \\ oc 
and 

Ti.± oc 
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dE 7 



(12) 



eV 



(E - Eft)(E - E*-) 



((E - E1) 2 + n 2 ){{E - El) 2 + n 2 ) 



dE. 



(13) 

The energy dependence of the real (red squares) and 
imaginary (black squares) parts of G < ^ for (E* << E^) 
and (E* w E') are shown in panels (al) and (a2) of 
Fig. [3 respectively. The corresponding bias dependence 
of Tj || (blue squares) and of T^j_ (green squares) are 
shown in panels (bl) and (b2) of Fig. [3 respectively. One 
can clearly see that the overall bias dependence of this 




FIG. 8: (Color Online) (a) Angular dependence of the net (a) 
field-like torque and (b) spin-transfer torque, for N c = 4 and 
various values of bias. 



simple model reproduces qualitatively that displayed in 
Figs[SJ: and[5Ji. More specifically, the switch-on behavior 
for Tj-ii and the sign change of Ti_± with bias are as- 
sociated with the relative position of the majority and 
minority QWS energies which lie within the bias energy 
window. When E* « E^, both components of the spin- 
transfer torque are dramatically enhanced, as shown in 
Fig. EIa2), due to the presence of higher-order poles in 
Eqs. (Ull) and (fl"3|). 

The net spin torque components 2]i(_n = J2i E c -^-11 (J-)' 
are not as strongly enhanced as the local torques, Tj ii 
and Ti which oscillate as a function of site i, due to 
the electron precession in the central FM (Fig[5]). The 
angular dependence of the net field- like torque, T±, and 
the spin-transfer torque Tji, is shown in Fig. [5^, and[5}3, 
respectively, for N c = 4 AS and for various bias values. 
For zero bias T± displays sinusoidal behavior, which how- 
ever, can change dramatically upon increasing the bias. 
The net field-like torque can be expressed^, [2J] as 



Ti 



-dE xc (9)/d9, 



(14) 
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FIG. 9: (Color Online) Bias dependence of the non- 
equilibrium effective bilinear, Ji, and biquadratic, J2, inter- 
layer exchange couplings for N c = 4. 



where 



E XC {0) = -Jicos{9) - J 2 cos 2 {B) 



(15) 



is the effective exchange coupling energy [251 ] . between 
and M^^). The out of equilibrium interlayer ex- 
change coupling has terms related to the spin-polarized 
tunnel current that can dominate and alter the coupling 
behavior under certain bias conditions}^. Here, J\ and 
J2 are the non equilibrium bilinear and biquadratic ef- 
fective exchange couplings, respectively, which are de- 
termined by fitting the angular dependence of T±(0) in 
Fig.[5]to the above expression for various biases. In Fig. [5] 
we show the bias dependence of J\ and J2 for N c — 4 AS. 
It is important to note that J\ (red circles) reverses its 
sign with bias, similar to Ti,±(V) in Fig.[6j:. On the other 



hand, the bias dependence of J2 (black squares) exhibits 
the switch on bias behavior found for Tjii(V"). Hence, 
there is a range of bias where Ji > Ju favoring perpen- 
dicular alignment of Mp and Mr.(g)[25fl. In contrast, the 
angular dependence of the net Xj | exhibits a skewed sinu- 
soidal behavior (J5J)) for any bias, similar to that found 
in spin valves. 



IV. CONCLUSION 

In summary, we demonstrate that the local spin- 
transfer torque can be dramatically enhanced, when the 
majority and minority QWS of different quantum num- 
ber are in close energetic proximity and lie within the 
bias energy window. This enhancement may in turn lead 
to reduction of the critical current necessary for CIMS in 
magnetic memories. The spin-torque enhancement crite- 
rion may be achieved by controlling the SPQWS through 
an external magnetic field or spin-dependent barriers. 
The SPQWS tune selectively the bias dependence of the 
spin-transfer and field-like components of the local and 
the net spin torque. This results to an anomalous angu- 
lar behavior of T± due to the bias interplay of the bilin- 
ear and biquadratic effective exchange couplings. Future 
work will be aimed to include the results for the local spin 
torques of these calculations as an input into the Landau- 
Lifshitz- Gilbert equation, to calculate the critical current 
for the CIMS. 
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